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Abstract 

We find new non-Abelian flux tube solutions in a model of Nf scalar fields in the fundamental 
representation of SU(N) x U(l) with iV < Nf (the "extended non-Abelian Higgs model"), and 
study their main properties. Among the solutions there are spinning strings as well as supercon- 
ducting ones. The solutions exist only in a non trivial domain of the parameter space defined 
by the ratio between the SU(N) and U(l) coupling constants, the scalar self-interaction cou- 
pling constants, the magnetic fluxes (Abelian as well as non-Abelian) and the "twist parameter" 
which is a non-trivial relative phase of the Higgs fields. 



1 Introduction 

Non-Abelian stringlike solutions have a long history which starts already in 1973 with the Nielsen- 
Olesen seminal paper pQ . Several general discussions were published [3 [H [3J 0] and explicit solutions 
(numerical of course) were obtained O [3 [H 02 [10l [11] for an SU (N) gauge theory with Higgs 
fields in the adjoint representation which completely break the symmetry. However, the major part 
of the activity in the field of cosmic strings |12j was concentrated on their Abelian counterparts. One 
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reason for this is that these non-Abelian string solutions have their flux directed in a fixed direction 
in the corresponding algebra so they arc essentially Abelian. 

In recent years, new kinds of non-Abelian strings were discovered during attempts to understand 
the phenomenon of confinement in QCD [T3l [T4J |T5l [16j HZ] , and their properties were studied 
extensively |18[ IT9|, [20] . These new solutions appear in models with a global (flavor) SU(Nt) 
symmetry in addition to the SU(N) local (color) symmetry based on scalar fields in the fundamental 
representation. They allow rotation of the non-Abelian flux in the Lie algebra, which makes them 
genuinely non-Abelian. When Nt > N these consist a generalization [2TJ[22] of the semilocal strings 
introduced originally within the extended Abelian Higgs model [22 [24] which is the Higgs system 
with a global SU(2) symmetry in addition to the local U(l). We therefore term the model discussed 
here the "extended non-Abelian Higgs model" . 

Most of the studies of these non-Abelian string solutions up to now have been limited to the self- 
dual (BPS) case. However, it is natural to go further and look for more general solutions as has just 
been done very recently [25[ 126] . This is the direction which we will take in this work, namely going 
beyond the BPS limit and it will be done together with allowing also for the possibility of rotation 
(i.e. spinning solutions) and of currents along the string axis. Spinning and superconducting cosmic 
strings have been found recently J37[ [2H] in the extended Abelian Higgs model which gives rise also 
to the (embedded) Nielsen-Olesen solutions. These new semilocal solutions, known as twisted, occur 
mainly outside the very peculiar (self-dual) limit of the coupling constants where the equations of 
the theory admit Bogomolnyi conditions. Another outstanding feature of twisted semilocal strings 
is that, when they exist, there is a continuous family of them, labelled by the "twist", a parameter 
entering through a space-time dependent relative phase of the Higgs field components. Note however 
the existence of the electro- weak superconducting strings |29] which exist without a "twist" . 

In the Abelian case, the local string is characterized by a magnetic field concentrated in a tube 
along the symmetry axis. Outside the core the magnetic field strength decays exponentially. The 
Higgs field vanishes on the axis and reaches asymptotically its symmetry-breaking value. 

For the twisted semilocal string the geometry is more involved. The magnetic and Higgs fields 
behave roughly as for local string but the configuration supports in addition an azimuthal ( "tangen- 
tial") component of the magnetic field. The source of this current is the additional Higgs component; 
its modulus is non-zero on the axis (forming a condensate) and vanishes outside the core and its 
phase is twisted. The effect of the non trivial phase can be appreciated once computing the gauge 
invariant Noether currents and the magnetic field. Figure [TJ gives a pictorial representation of the 
effect of the twist. 

In this work we will present the non-Abelian analogues of the Abelian semilocal twisted strings 
and study their properties like current and angular momentum. Untwisted purely magnetic local 
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Figure 1: The magnetic field and few spiral field lines of a twisted Abelian semilocal string. The thick line 
is the string axis. 

non-Abelian strings will be discussed briefly as a special case. The field equations corresponding 
to the models presented here are non linear and coupled and do not admit explicit solutions. We 
therefore rely on numerical techniques to construct the solutions and calculate the physical quantities. 
Several figures are necessary to illustrate the extreme richness of the solutions. 

The existence of these new kinds of strings raises the question of the nature of the gravitational 
fields of these strings and the possibility of new features in this respect. Some initial work has been 
already done in this area [30j - still for the BPS case only and we will turn to that question in a 
future publication. Another issue which is beyond the scope of this work is the effect of spin on the 
reconnection probability [311 E2] of non-Abelian strings. 

This paper has the following plan. In section 2 we present the extended non-Abelian Higgs model 
and discuss its relation with its Abelian counterpart. In section 3 we derive the field equations 
and obtain the physical quantities which are used to characterize the solutions: energy angular 
momentum, currents and charges. In section 4 we present the various string solutions and discuss 
their properties across their parameter space. Section 5 contains our conclusions. 

2 The extended non-Abelian Higgs model 

Our general framework is based on a Lagrangian describing a multiplet of Nf scalar fields with local 
invariance under SU(N) x U(l) and a global invariance under SU(Nf). The local symmetry further 
requires iV 2 — 1 non-Abelian gauge fields and one Abelian field. The scalar fields can be written in 
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term of a matrix with elements $ as where 1 < a < N and 1 < s < Nt transforming according to 

& br = U ba (x)$ as Sl (2.1) 

where U and S are matrices in the fundamental representations of SU(N) and SU(Nf) respectively. 
The Lagrangian is: 

C = (D^ as )*(D^ as ) - V($ as ) - \F; v F a ^ v - \f^F^ (2.2) 
The standard definitions are used for the covariant derivative and gauge field strengths : 

D^ s = (d„ - ie x A^ - ie 2 A«T a ^ $ s (2.3) 

F va , = d vt A v -d v A ll ; F% v = d„Al - d v A« + 2e 2 f ahc A\Al (2.4) 

where f abc are the structure constants of the gauge group and ■^T a are the (hermitian) generators 
in the fundamental representation. We use the normalization Tr(T a r b ) = 28 ab . 
The most general renormalizable symmetry breaking potential is 

Be 2 

V = ae?(Tr($ t $) - Nv 2 ) 2 + -J-Tr(&T a <£>)Tr(&T a <S>) (2.5) 

where a , B and v are positive constants. It can be written also as: 

V = ae?(Tr($ 1 '$) - Nv 2 ) 2 + Bel \ Tr(&$ $ f $) - -^(Tr^®)) 2 ^ (2.6) 

The first term forces the scalar fields to develop a non-trivial minimum, while the second term forces 
the minimal configurations <f> to be such that $T$ = v 2 j 

This general system contains some well-known special cases, or seen from the other direction, 
may be considered as a generalization of several systems. From the point of view of the present 
discussion our model generalizes the extended Abelian Higgs model which corresponds to N = 1 
(and e2 = 0, while Nf is arbitrary) in our terminology. In the Abelian case, the second term of 
the potential vanishes identically. This model allows for two kinds of stable string solutions as was 
first discovered [53] for the case Nf = 2. The first is just the embedded Abelian (Nielsen-Olesen) 
flux tube which in these circumstances is stable [33J for a < 1/2 . The second is sometimes called 
"skyrmion" because of the relation with the er-model lumps [33] and exists only in the self-dual 
limit where the masses of the Higgs and the gauge particles are equal (or a = 1/2). Its stability 
is guaranteed by a Bogomolnyi-type argument. Both solutions are termed semilocal strings and 
their discovery was a surprise at the time, since the vacuum manifold of the model is a (2Nf — 1)- 
dimensional sphere which does not give rise to non-contractible loops. This discovery gave an explicit 
example for a system where the non-triviality of the first homotopy group of the vacuum manifold is 
not a necessary condition for the existence of stable string-like solutions. In this system, the kinetic 
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(gradient) term plays a crucial role since it has a lower symmetry than the potential term, and it 
is this "mismatch" between the different symmetries which enables the existence of these solutions 

[33 [Ml \M- 

The two kinds of solutions mentioned above are static, but stationary spinning semilocal strings 
also exist in this system, as well as solutions which carry a persistent current. These solutions 
[2"Tl |2"8] belong to a new family of solutions named twisted semilocal strings. The "twist" is realized 
by a relative phase between the two Higgs field components, which changes along the string axis 
and may be also time dependent. These twisted strings exist for a > 1/2 and form a continuous 
family parametrized by this "twist" parameter. The main physical property of these twisted strings 
is the current along the string axis which exists in both the static and stationary cases and gives 
rise to azimuthal component of the magnetic field. Their energy in their rest frame decreases with 
growing current which implies that they are energetically more favorable to be produced during 
(cosmological) phase transitions than the embedded Abelian strings. These new solutions bifurcate 
with the embedded Abelian flux tubes in the limit of vanishing current, thus clarifying the nature 
of the "magnetic spreading" instability [33, 35] of the embedded Abelian flux tubes for a > 1/2. 



3 Non- Abelian Semilocal Strings 



3.1 Field Equations for Stringlike Solutions 

We look for cylindrically symmetric non-Abelian stringlike solutions and for definiteness, we limit 
ourselves to the case Nf = N + 1. But, the ansatz presented below and the corresponding equations 
can be generalized easily. 

The electromagnetic potential and its non-Abelian counterpart will be: 

A^dx" = A (r)dt + A g (r)d9 + A z (r)dz ; r° A°dx>* = (l (r)dt + A e {r)dO + A z {r)dz^j G (3.1) 

where G belongs to the Cartan subalgebra of the gauge group. For SU (2) there is only one possibility, 
i.e. G — (T3 but already for SU(3) which is a rank 2 group there is more freedom and G is a linear 
combination of the two (commuting) diagonal Gell-Mann matrices. We will denote the diagonal 
elements of G by G a with a = 1, N. 
For the scalar field we take: 

/ <^(r)e iai •••0 \ 

if 2 {r)e ia2 ••• 



$ = 



(3.2) 



ip N (r)e laN , p N (r)e iXN J 
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which follows (and actually generalizes) the form of [2] or [3D]. The phases a, and xn will have 
linear dependence on t, 9 and z: an = u)(t + rriiO + kiZ and xn = t&t + n6 + qz. 

To write the field equations, it is convenient to use a dimensionless coordinate, to rescale the 
fields and "phase parameters" appropriately and to define a relative gauge coupling constant: 
x , , , , e 2 



wei ei 
(A,(r),A,(r))=«(>4o(a:),A,(a;)) , (A (r), A z (r)) = —(M^), A,(x)) 

&2 



(3.3) 



With these dimensionless quantities, the field equations of the system are written below. For the 
scalar fields we obtain (recall: G a are the diagonal elements of G) : 

(x<P'a,y ( ( m a - M- AgGg) 2 ~ 2 fo A\ 

— - — =| — (uj a - A - A G a ) + (k a - A z - A z G a ) \ fa (3.4) 

ip a , a < N 



+2 



V 6 / \ b^a 



(xip' N y _ j (m N -Aq- A e G N f 



+2 



(lu n -Ac- AqG n ) 2 + (fejv - A z - A Z G N ) 2 J ^ ( :■{.:. ) 



AT-l 



b=l 



{xp' N )' ( (n- A e - A e G N ) 2 



(zu-A - AoG N ) 2 + {q-A z ~ A Z G N ) 2 p N (3.6) 



+2 



where there is no summation on repeated indices and all sums are for a, b = 1, N unless indicated 
otherwise. For the Abelian-gauge fields, we have 
(xA'J 



-2 EK - A - A G a )^ - 2(ro - A - A Gat)p^ 

a 

z{A' e /x)' = -2 ^(m a - Ag - A e G a )<p 2 a - 2(n - A 9 - A e G w )p^ 

a 

= -2 — A z — A,G a )^ - 2( q -A z - A Z G N ) P 2 N 



Finally, for the non-Abclian fields 

(xA' y 



= -2S 2 I J2^ a ~ A °- A Ga)<p 2 a G a + {w - A - A G n )p 2 n Gn 

\ a / 

iA' e /x)' = -2S 2 (E( m « - A e - A*G a )<plG a + (n-A e ~ A e G N )p 2 N G N 



(xA' z y 



= -25 2 J2( k - ~ A *~ AzGa^lGa + (q-A z - A z G N )p 2 N G N 



(3.7) 
(3.8) 

(3.9) 

(3.10) 
(3.11) 
(3.12) 



In order to obtain localized solutions, appropriate boundary conditions should be imposed. The 
boundary conditions will guarantee that the solutions will be regular on the axis x — and approach 
a vacuum for x — > oo. Regularity at x = yields 

p«(0) = , A g (0) = , A e (0) = , /^(O) = , A^(Q) = , 1^(0) = (3.13) 

while the requirement for asymptotic approach to the vacuum is somewhat more involved. For the 
scalar field we simply have 

^ a (oo) = 1 , pat(oo) = (3.14) 

which makes a clear distinction between the N fields tp a and the pjy field. This last one will either 
trivially vanish, or have a unique behavior of starting at a non- vanishing central value and decreasing 
monotonically to zero. Therefore, to avoid singularity at x — all the solutions discussed in the 
following will have n = 0. 

For the asymptotic behavior of the gauge fields we start by the observation that under a Lorentz 
boost in the z direction, the quantities (u>i,ki) transform as Lorentz two- vectors. Since we will 
consider in this paper only the case when these two vectors are spacelike (corresponding to the 
magnetic case in the classification of [28]), we can choose one of these two- vectors to be of the form 
(0, ki) without loosing generality. Independently, one can perform residual gauge transformations of 
the form 

U res = expi((at + bz)I + (aft + b'z)G) (3.15) 

and assume in the following 

A (oo) = , A z (oo) = , i (oo) = , A z (oo) = 0. (3.16) 

As for the azimuthal components of the gauge fields we have to impose 

m a ~A e {oo) -A e (oo)G a = (3.17) 

These are iV equations for the two unknowns Ag(oo) and Aq(oo), so the three magnetic numbers m a 
cannot be independent. Moreover, since they must be integers, not every G in the Cartan subalgebra 
is possible. A closed formula for these relations is rather complicated so we demonstrate this by two 
examples for the case N — 3 where we can parametrize G by: G — A3 cos?/> + Ag sin?/; where A, are 
Gell-Mann matrices. For tp = tt/6 we have the relation mi = and A$(od) = (mi —ms)/y3, while 
for -0 = 7r/3 we have 7712 = (mi + m^)/2 and accordingly Ag(oo) = (mi — m3)/2. Ag(oo) however is 
given always by ^4#(oo) = (mi + rri2 + m^,)/Z. The concrete solutions which we will present will be 
mainly with N = 3 (other values of N can be treated similarly) . 
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3.2 Physical quantities: energy, angular momentum, currents and 
charges 

In this section, we express several quantities which characterize the different types of solutions. 
The quantities are presented for string solutions of the extended non-Abelian Higgs model. The 
corresponding formulas for the "non extended" model (i.e. Nf — N) are obtained by setting pn = 
and Xn = in Eq. (|3.2p and those that follow. 

The energy and angular momentum (per unit length) of the solutions are obtained by computing 
the integrals 

POO />OC 

E = 2ttv 2 / dxx(B +£i+V) , J = 2ttv 2 / dxxj (3.18) 
Jo Jo 

The energy densities £q,£i,V, correspond respectively to the scalar part, the gauge field part and 
the potential. They are given by 

, 2 . / (m tt - Ag- A g G a ) 2 



9a 



+ (w -Aq- AaGaY + (k a -A,- A z G a ) 2 <pi 



(3.19) 



-Pn 



, 2 / (n-A e - A e G N f 



+ (u7-A - A Gn) 2 + (q-A z - A Z G N ) 2 p 2 N 



B t = \ {A' 2 + (A' e /x) 2 + A' 2 ) + ± (A» + (A'/x) 2 + A?) 



(3.20) 



V a 



+ P 2 N -N\ +(3 



a \ a / 



The angular momentum density J has the form 

J = 2j2(^a - A - A G a ){m a - A e - A e G a )ip 2 a 



(3.21) 



(3.22) 



+ 2{w -Aq- A G N )(n - Ag - A e G N )p 2 N + A' A' g + —A' A'g 

The solutions are also characterized by the Noether charges and currents associated with the 
global SU(N + 1) symmetry of the Lagrangian. The charges (per unit length) are the integrals of 
the time components of the current densities 



MR) = -* {KpRpqD^aq - (D^ ap )*R pq ^ aq ) 



(3.23) 



where R denotes any of the generators of the global symmetry and the summation convention is 
used for both the local and global indices. For the "diagonal" ansatz, Eq. (|3.2II . we get the following 
expressions for the global charge densities 



Jo 



{R) = 2J2 Ra(uj a -A Q - A G a )^ 2 a + 2R N+1 {vo - A - A G N )p 



(3.24) 



where we denote by R s the diagonal elements of the commuting flavor generators as we do for local 
symmetry. Similarly for the current densities: 

j z {R) = 2^R a {k a -A z -A z G a )<pl + 2R N+1 {q-A z -A z G N )p 2 N . (3.25) 

a 

Accordingly, the charges and currents will be 

,.00 JV+1 .oo JV+1 

Q(R) = 2tt / dxxjo(R) = R *Q S > *( R ) = 277 / dx x] z {R) = RsIs ( 3 - 26 ) 

where the "single field contributions" Q s and I s are defined in an obvious way. In terms of these, 
the angular momentum J of the string can be written as J = ^ m a Q a + tiQn+i- 

Since the non-Abelian charges all vanish due to the boundary conditions, these "single field con- 
tributions" are not independent but satisfy the following identities among the charges and currents 

JV+1 JV+1 N N 

J2Qs=J2 I s = ° ' G -Q- + G JvQjv+i = G " Ia + G nIn+i = 0. (3.27) 

s—1 s—1 a—1 a—1 

These identities allow one to express two of the currents (or charges) in terms of the others. As a 
consequence, for an SU(N + 1) global symmetry, the N diagonal global currents can be expressed 
in terms of N — 1 independent ones. For example, for local SU(2) and global SU(3), these relations 
lead to 1\ = 0; therefore I2 — —I3 by the left equation of (|3.27p and the two diagonal global currents 
are proportional to I3. 

For local SU(3) and the choice ip = 7r/6, the global SU(A) currents depend on two of the "single 
field" currents / s say, I2 and I4. Actually two of them are just proportional to I2 and I4 and the 
third is a combination: 

I{R 3 ) = -2I 2 ; I{R & ) = -±={3h + 2h) ; I(R l5 ) = -\J\h- (3.28) 

These identities provide useful crosschecks of the numerical solutions discussed next. 



4 Discussion of the solutions 

We solved the system of non- linear equations (|3.4[) - (|3.12[) by using a numerical routine based on the 
collocation method 35J. The solutions were constructed with grids involving typically 1000 points 
and with accuracy of the order of 10 -8 . The solutions we found are described in the following. 

4.1 Special Case: Local strings 

Before discussing the semilocal strings, we pause to describe the simplest kind of non-Abelian ones, 
namely the local non-Abelian strings. They are obtained by truncating the general system to the 
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Figure 2: A non-Abelian SU(3) local string: Profiles of the three Higgs fields tpi, tpz, ip3, the two gauge 
fields Ag and Ag (A' g in the figure) and the corresponding magnetic fields for mi = 2, m2 = 1, mz — 1 . 
This choice yields the relation tp2 = y>3- The Abelian flux is 3/2; the non-Abelian flux is 1/V3. The relative 
gauge strength is S — 1 and the parameters in the potential are: a = /3 = 1. 

case with Nf — N so pn = and \n = 0. This requires the further substitutions za = n = q = in 
the field equations above. 

Inspection of the boundary conditions for the gauge fields as x — ► oo, reveals that in the gauge 
choice p.!6p all the parameters u)j, ki should vanish, leading to Aq(x) = A z (x) = Aq(x) = A z (x) = 0. 
So we solve equations (|3.4l) - (|3.12p with the above substitutions and with the following boundary 
conditions : 

p o (0) = 0, Ae(0) = 0, A e (0) = (4.1) 
^ O (oo) = l, m a - A e (oo) - Ae(oo)G a = (4.2) 

The only stringlike solutions in this system are of the magnetic type, labelled by the integers m a 
which are consistent with the value we pick for the angle ip in the 3-8 plane of the Cartan subalgebra. 

In Fig. [5] we show the field profiles of a typical SU(3) local string solutions with tp = tt/6. The 
Abelian flux is still quantized but is not an integer as can be easily seen from the figure. Its value 
is 3/2. Similarly, the non-Abelian flux is 1/V3- 
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4.2 Semilocal strings 

Now we return to the general case of the full system with Nf = N + 1, that is equations (|3.4|) - (|3.12|) 
with non- vanishing pjj and \n- The solutions have in general (zu, q) ^ while we can still exploit 
the remaining symmetry to require again ujj = kj = for j = 1, . . . , N. 

In contrast with the local case, the equations for the fields A , A z , A and A z are not trivially sat- 
isfied if pn(x) ^ and we have now non-trivial boundary conditions for these gauge components. The 
boundary conditions are contained in equations (|3. 13|) . (|3.14p . (|3 . 1 6|1 and (|3.17j) and we will not re- 
peat them here. We just point out that the scalar and gauge fields reach their asymptotic values with 
exponential corrections. In particular, the field pjv decays according to pjv oc exp [— \J q 2 — zu 2 x], 
which explains why this field can decay slower than the others for \J q 2 — w 2 < < 1 . We have there- 
fore to solve for p^(x), Aq(x), A z (x), Ao(x) and A z (x) in addition to the N scalar components and 
the two gauge components Ae and Ag. 

Along with [57j, we call solutions of this type "twisted" solutions. The parameter q will be referred 
to as the "twist" parameter. The solutions will be generally "twisted", but special "untwisted" 
solutions (with q = 0) also occur. We also distinguish between static solutions with zu = and 
stationary for zu ^ 0. Their properties will be studied below. 

4.3 Purely magnetic ("untwisted") solutions 

The simplest kind of semilocal solutions is the embedded (non-Abelian) local flux tubes. These are 
purely magnetic solutions of the full system with pjv = 0. It is therefore evident that A$ = A z = 
Aq = A z = and the non- vanishing fields, the N scalars ip a , and the two gauge components Ag and 
Ag, are the same as for the local strings of Sec. 14.11 

Another kind of purely magnetic strings is the untwisted semilocal strings which generalize the 
"Skyrmion" solutions of the Abelian model mentioned in Sec. [2] These solutions are only self-dual 
as in the Abelian case, and are contained already among the results of earlier studies [211 E21 SO] . 
In our parametrization they exist for a = 1/2, /3 = 1, 5 = 1. 

The more general kinds of solutions which contain "twists" [27l [28] will be discussed in the next 
sections. 

4.4 Twisted static and stationary Semilocal strings 

The space of twisted semilocal strings decomposes into three regions which are labelled according to 
the sign of the norm of the two- vector (zu, q), that is according to the sign of zu 2 — q 2 . In the "timelike" 
region (zu 2 — q 2 > 0) we were unable to find localized solutions as happens also in the Abelian case 
[28j (recall also the asymptotic behavior mentioned above), so we do not discuss them further. 
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Solutions with vo 1 — q 2 = are known as "chiral" and will be discussed below among the stationary 
solutions. If we assume the two- vector (zu, q) to be spacelike, we can set w — by an appropriate 
Lorentz boost. The field equations for Aq and Aq are satisfied trivially by Aq(x) — Aq(x) — and 
all the charges Q s vanish with this choice. The components A z and A z as well as the currents I s do 
not vanish. 

In the Abelian case (N = 1 and Nf = 2 with irrelevant (3) [S7J 12H] twisted semilocal strings 
appear for large enough values of the potential coupling constant a (precisely for a > 1/2 in our 
notations). They form continuous family of static solutions and are labelled by the "twist parameter" 
q appearing in the phase of the scalar function p^\ this parameter takes values in a finite interval, 
that is q € [0, q cr ] where q cr is the maximal twist where bifurcation with the embedded flux tube 
occurs. This value of q cr depends on the coupling constant a. In the limit a — > the function pi 
converges uniformly to the null function and the branch of twisted solutions bifurcates into the local 
string solutions. 

Here, we consider the extended non- Abelian Higgs model with our potential (|2.5[) and find that 
analogous string solutions exists when the potential parameters, a, (3, are chosen large enough. A 
similar critical phenomenon occurs, and in addition to the critical twist there exists now a critical S. 

First, we study the case of an SU(2) gauge group and a global SU(3). We find twisted semilocal 
strings corresponding to different choices of the magnetic parameters (mi, m-i). Fig. [3] shows the 
field profiles of two typical static twisted semilocal strings with different windings and therefore 
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Figure 3: A non- Abelian SU(2) static twisted semilocal strings: Profiles of the three Higgs fields <pi, </?2, pi 
and the components of the two gauge fields Ag,A z and Ag, A z (A' g , A' z in the figure). The other parameters 
are: a = 1, (3 — |, 8 — 1. Left: Note that the non-Abelian magnetic flux vanishes and that A z = — A z . The 
Abelian flux is 1. Right: Note that A z = —A z . The Abelian flux is 3/2; the non-Abelian flux is 1/2. 
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Figure 4: A non-Abelian SZ7(3) static twisted semilocal string {A' stands for A). 5 = 1. 

different fluxes. Note the additional scalar field pi which decreases slowly, yet still exponentially. 
The relations among the currents, Eq. (|3.27[) yield in this case 1\ = (see explanation below (|3.27p ) 
and we observe that it is consistent with the numerical results giving A z = —A z for this case as is 
indeed seen in Fig. [3J 

Analogous twisted solutions were obtained for N = 3 and Nf — 4. Fig. 0] depicts a typical 
solution corresponding to the specific Cartan generator fixed by tj) = tt/6. The magnetic fields 
along the string axis (Abelian and non Abelian) and the energy density of a twisted solution (with 
l°3(0) = 1.25 or q = 0.34) are compared in Fig. [5] with those of an embedded local string. Note 
that the twisted strings have a non-vanishing tangential magnetic field (Bg-see figure [1} whose 
contribution to the flux vanishes. 

The solutions constructed above are static but they can be made stationary by applying a boost 
along the string axis (z). Alternatively, stationary solutions can be obtained directly by setting the 
parameter w non zero in the field equations. Interestingly, allowing the parameter tn 7- confers the 
strings with angular momentum as was first demonstrated in the Abelian case |27[ I28| . The resulting 
angular momentum is therefore more like the "orbital" kind rather than an intrinsic spin-like. 

The profiles of twisted solutions corresponding to mi — 2, m,2 — TO3 = 1 and q — 0.5 are presented 
in Fig. [6] for w = and w — 0.495 (i.e. very close to the chiral limit vj = q). We note that the 
chiral limit vj — > q leads to a different kind of solutions with p 3 decaying as a power instead of 
exponentially. This tendency is apparent in Fig. [S] namely, through the decay of the function p%. 
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Figure 5: The energy density and magnetic fields along the string axis (B' stands for B) are superposed for 
a twisted semilocal string P3(0) = 1.25 or q — 0.34 (solid lines) and an embedded local one (dashed lines). 




Figure 6: Field profiles (A' stands for A) of two N = 3 twisted solutions with a = /3 = S = 1 , rh = (2, 1, 1), 
q = 0.5. Solid lines: static (w = 0); dashed lines: stationary with w = 0.495. 

4.5 Analysis of parameter space 

Since we know from the Abelian case (N = 1) that twisted semilocal strings exist only is a restricted 
domain of the parameter space, it is natural to try to figure out this domain in the present non- 
Abelian model. The Lagrangian (|2.2[) is characterized by four parameters: a, (3 of the scalar potential 
and two gauge coupling constants e\ and e%. However, only the relative strength 5 = e-xje\ appears 
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in the field equations so the parameter space of the model is three-dimensional. 





Figure 8: Energy, global currents (see (|3.28|0 . central magnetic fields (B' stands for B) and central "twisted" 
scalar field as a function of the coupling ratio S 2 for static twisted semilocal strings with q — 0.4 (solid lines) 
and embedded local ones (dashed lines). Both for a = f3 = 1 and in = (2, 1, 1). Note that the twisted string 
has lower energy. 
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The parameter space of the solutions is far more complicated since it contains also the twist q, 
the other phase parameter w and the flux numbers. 

We investigated this problem in the case of an SU (3) gauge group and we present here the main 
results in various planes in this higher dimensional parameter space. We took representative values 
of a = = 1 for which it turns out that twisted semilocal strings exist. 

The effect of changing the parameter 6 on the twisted solutions is illustrated in Fig. [7] where the 
critical g cr -value is plotted as a function of S 2 for two different choices of to = (mi, TOa, TO.3). The 
choice (2, 1, 1) corresponds to ip = t/6 and the choice (3, 2, 1) to ip = ir/3 as explained at the end of 
Sec. 13.11 The twisted semilocal strings exist only for q < q cr and for values of S 2 below a maximal 
value of S 2 r . This value depends on to and determines a domain of existence of twisted solutions in 
the (S, q) plane. 

In Fig. [H a few physical quantities characterizing the solutions corresponding to q = 0.4 are 
presented. We find 5 2 r = 5.55 which is consistent with the value that is obtained from Fig. [7] We 
further studied the critical phenomenon appearing when the maximal value of d (with q fixed) is 
reached as seen in Fig. [5J Our numerical results suggest in particular that the twisted solution 
bifurcates into the embedded local string at S = 8 cr . Note especially the vanishing of I 2 , I4 and 
p 3 (0) as 5 — > S cr . 

The evolution of the physical parameters characterizing the twisted solutions for 8 — 1 is further 
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Figure 10: Evolution of the Mass and the magnetic fields (B' stands for B) on the axis as functions of zu 

presented in Fig. [9] as a function of the twist parameter q. We see from the plot that the solution 
bifurcates into the embedded local string for q — 0.597 which is of course consistent with q cr (l) 
of Fig. [7J We notice that the central value of the non Abelian magnetic field, -B(O) changes sign 
between the bifurcation point and the low twist region (i.e. q — > 0). 

The effect of a boost along the string axis, which is reflected by a non- vanishing w, on some 
physical quantities is presented in Fig. \W[ In this figure E, J, B(0) and B(0) are plotted versus 
w for a fixed value of q (for instance q = 0.5) and for the two values i/j = tt/6 ,7t/3 which require 
different sets of m a as explained above. We stress that the plots are made for fixed q, so the solutions 
involved are not related to each other by boosts along the z direction. In particular, the magnetic 
fields on the z-axis vary with w. Note also that the chiral limit has nontrivial consequences on the 
mass and spin of the solutions (per unit length) which diverge in this limit as indicated in Fig. [TDJ 



5 Conclusions 

This work presents string solutions in the extended non- Abelian Higgs model which may be seen as 
an expansion of recent discussions of Abelian twisted semilocal strings [S7J HH] and of non- Abelian 
untwisted ones [2TJ [22l [25l [26] . The solutions are classified according to the twist parameter q as 
twisted or untwisted (for q = 0). The twisted solutions are characterized by a persistent current 
along the string axis. The twisted solutions are further classified as static {vj = 0) or stationary, 
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where the latter are also characterized by a non-zero angular momentum. The extended non-Abelian 
Higgs model allows also two kinds of special solutions - both untwisted and both purely magnetic: 
The embedded local non-Abelian flux tubes and the self-dual "Skyrmions" . 

The occurence of Twisted semilocal string and their corresponding conserved currents is closely 
connected to the existence of a scalar field which condensates on the axis and vanishes asymptotically. 
For Nf = N all scalars are Higgs-like fields which acquire asymptotically a non vanishing expectation 
value due to the potential. As a consequence, there is no place for a condensate and no twisted 
semilocal string are possible in this case. 

All those solutions were constructed by numerically solving the field equations with the appro- 
priate boundary conditions. The solutions exist only on a non trivial domain of the parameter space 
of the possible solutions. We have tried to determine this domain qualitatively and obtained a lot 
of information about the pattern of the solutions. We demonstrated the bifurcation phenomenon 
which occurs in this case at a critical value of the twist q which now depends on the various fluxes 
and on 8 = in addition to the parameters in the potential. 

A crucial step for deciding the physical relevance of non-Abelian twisted semilocal strings would 
be to study their stability by performing a systematic normal mode analysis, as was done recently 
[36j for the extended Abelian Higgs model. Another important issue is the understanding of the 
low-energy dynamics of the non-Abelian twisted semilocal strings. This can be done e.g. by using 
the geodesic approximation of |37[ 138] which is based on a parametrization of the zero modes in 
terms of collective coordinates. For untwisted local and semilocal strings, this analysis is reported in 
details by Aldrovandi (see sec. VI of [30]). Since the phases of twisted semilocal strings depend on t 
and z, the parametrization of the orientational moduli in terms of slowly varying functions of these 
coordinates does not seem to be straightforward. The parametrization of the fluctuations used in 
[36j might turn out to be very useful. We plan to address these questions in a further publication, 
as well as the gravitating counterparts of twisted strings. 

Another direction is motivated by the non-Abelian superconducting strings, recently obtained 
by Volkov [29j within the SU(2) x U(l) electroweak theory. These "untwisted" strings can be used 
as a starting point towards constructing twisted solutions which will be related to a non-Abelian 
subgroup of the gauge group. 
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